Coherent frequency conversion in a superconducting artificial atom with two internal 

degrees of freedom 



F. Lecocq, I. M. Pop, I. Matei, E. Dumur, A. K. Feofanov, C. Naud, W. Guichard, 0. Buisson 

Institut Neel, C.N.R.S.- Universite Joseph Fourier, BP 166, 38042 Grenoble- cedex 9, France 

(Dated: January 20, 2012) 

By adding a large inductance in a dc-SQUID phase qubit loop, one decouples the junctions' 
dynamics and creates a superconducting artificial atom with two internal degrees of freedom. In 
addition to the usual symmetric plasma mode (s- mode) which gives rise to the phase qubit, an 
anti-symmetric mode (a- mode) appears. These two modes can be described by two anharmonic 
oscillators with eigenstates \n s ) and \n a ) for the s and a- mode, respectively. We show that a strong 
nonlinear coupling between the modes leads to a large energy splitting between states |0 S , l a ) and 
|2 s ,0 a ). Finally, coherent frequency conversion is observed via free oscillations between the states 
|0 a ,l o ) and |2 s ,0 a >. 

PACS numbers: 



In atomic physics the presence of multiple degrees of 
freedom (DoF) constitutes a precious resource for the 
development of quantum mechanics experiments. In 
trapped ions or nitrogen-vacancy centers, the V-shaped 
energy spectrum enables very high fidelity readout of the 
states encoded in the first two levels by using the fluores- 
cence properties of the transition to the third level [l], Q . 
In quantum optics, the multiple DoF of the optically 
active crystals have lead to many quantum effects such 
as Coherent Population Trapping (CPT) and the associ- 
ated effect of Electromagnetically Induced Transparency 
(EIT) 0, spontaneous emission cancelation via quantum 
interference 0, [H[ or generation of entangled photon pairs 
0| . In the field of superconducting qubit s, experimental 
efforts have mainly focused on two- level systems @,[s[ and 
multilevel systems with a single DoF. Supercon- 

ducting artificial atoms currently need additionnal DoF 
in order to realize A, V, N or diamond-shaped energy lev- 
els and therefore to perform new quantum experiments 
[3, [Hj. Only recently a superconducting device with 
y-shape energy spectrum was experimentally considered 
[l6| and a three DoF superconducting ring was developed 
for parametric amplification [17 ]. 

In this letter we present an artificial atom with two 
internal DoF which constitutes a basic block for the re- 
alization of V- or Diamond-shaped energy levels. In this 
system we benefit from the natural nonlinear coupling be- 
tween the two DoF to observe a coherent frequency con- 
version process in the time domain. Contrary to previous 
frequency conversion proposals and implementations in 
solid state devices, we observe this process in the strong 
coupling limit, where multiples oscillations can be seen 
before losing coherence [18|, LL9| , and without any exter- 
nal coupling device or additionnal source of power (2o| . 
In addition this system could be used for triggered and 
high-efficiency generation of entangled pairs of photons, 
a key component for quantum information (2lj . 



The circuit is a camelback phase qubit [22| with a large 
loop inductance, ie a dc-SQUID build by a supercon- 



ducting loop of large inductance L interrupted by two 
identical Josephson junctions with critical current I c and 
capacitance C, operated at zero current bias (see Fig|T|. 
As we will see in the following, the presence of a large loop 
inductance modifies dractically the quantum dynamics 
of this system. The two phase differences <\>\ and <p2 
across the two junctions correspond to the two degrees 
of freedom of this circuit, which lead to two oscillating 
modes: the symmetric (s-) and the anti-symmetric (a-) 
plasma modes [23[ . The s- mode corresponds to the well- 
known in-phase plasma oscillation of the two junctions 
with the average phase x s = + 02)/2, oscillating at 
a characteristic frequency given by the plasma frequency 
of the dc-SQUID, uo^. The a- mode is an opposite-phase 
plasma mode related to oscillations of the phase differ- 
ence x a = (0i — fa) /2, producing circulating current 
oscillations at frequency uu^. In previous experiments 
[sl EH Hll , the loop inductance L was small compared 
to the Josephson inductance Lj = ^o/(2ttI c ). There- 
fore the two junctions were strongly coupled and the 
dynamics of the phase difference x a was neglected and 
fixed by the applied flux. The quantum behavior of 
the circuit was described by the s-mode only, showing 
a one-dimensional motion of the average phase x s . Here- 
after we will consider a circuit with a large inductance 
(L > Lj) that decouples the phase dynamics of the two 
junctions. This large inductance lowers the frequency of 
the a-mode and the dynamics of the system becomes fully 
two-dimensional. The a-mode was previously introduced 
to discuss the thermal and quantum escape of a current- 
biased dc-SQUID 24], 25| but its dynamics was never ob- 
served. We present measurements of the full spectrum of 
this artificial atom, independent coherent control of both 
modes and finally we exploit the strong nonlinear cou- 
pling between the two DoF to observe a time resolved up 
and down frequency conversion of the system excitations. 

The electronic properties of a symmetric and inductive 
dc-SQUID at zero-current bias can be described by a fic- 
titious particle of mass m = 2C(&o/2tt) 2 moving inside 
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FIG. 1: Description of the device, (a) A micrograph of 
the aluminium circuit. The two small squares are the two 
Josephson junctions (enlarged in the top right inset, 10/ira 2 
area, I c = 713 nA and C = 510 £F) decoupled by a large in- 
ductive loop (L = 629 pH). The width of the two SQUID arms 
were adjusted to reduce the inductance asymmetry to about 
10%. Very narrow current bias lines, with a large 15 nH- 
inductance, isolate the quantum circuit from the dissipative 
environment at high frequencies . The symmetric and anti- 
symmetric oscillation modes are illustrated by blue and red ar- 
rows, respectively, (b) and (c) Potentials of the s and a-mode 
respectively, for the bias working point (J& = 0,$b = 0.374>o). 
(d) Schematic energy level diagram indexed by the quan- 
tum excitation number of the two modes \n 3 ,n a ) at the same 
working point. Climbing each vertical ladder one increases 
the excitation number of the s-mode, keeping the excitation 
number of the a-mode constant. The two first levels, |0 s ,0 a ) 
and |l s ,0 a ), realize a camelback phase qubit. 
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where Ej = $>qI c /(27t). Hereafter we will consider the 
particle trapped in a local minimum (x™ in ,x™ in ) veri- 
fying xf in = and x™ in + (L/2Lj) smx™ in = tt<V$ . 
By expanding the potential at this minimum up to fourth 
order, its quantum dynamics is given by the Hamil- 
tonian H = % s + % a + £s,a where % s and % a are 
anharmonic oscillator Hamiltonians describing respec- 



tively the symmetric and antisymmetric plasma modes 
[26L 27]. C S}a describes the coupling between the two 
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oscillators. H a = huj^ \{p 2 a 
where a = s, a with (ojp) 2 = (2E j / m) cos x 7 ^ 
(uj%) 2 = (uj s p ) 2 + (4Ej/m)(Lj/L). The operators x a and 
p a are the reduced position and momentum operators in 
both directions. The coupling term has a very simple 
expression at zero current bias: 



C s ,a = hg 2 ix 2 s x a + hg 2 2X z s x z a 



C.2 -2 



(2) 



with %2i = —Ej(h/m) 2 (cj^y/oJa) sinx^ m and % 2 2 = 
-Ej{fi/m) 2 {uj s v uj^)- 1 cosx™ n /2. In the following, we 
will define \n 3 ,n a ) = \n s )\n a ) as the eigenstates of uncou- 
pled hamiltonian H s -\-H a ^ where \n a ) indexes the energy 
levels of each mode. The potentials associated with the 
5-mode and a-mode are depicted in the figure FigJTjD and 
FigHt for a working point (I b = 0, $ b = 0.37$ ). At the 
same bias point the complete spectrum of the system is 
presented in Fig JUL 

The complete aluminum device is fabricated using an 
angle evaporation technique without suspended bridges 
[281 ]. and it is presented in FigJTJi. The measurements 
were conducted in a dilution refrigerator at 40 mK us- 
ing a standard experimental configuration, previously de- 
scribed in Ref.Q. 

The readout of the circuit is performed using switch- 
ing current techniques. For spectroscopy measurements 
we apply a microwave pulse field, through the current 
bias line (see Fig , followed by the readout nanosec- 
ond flux pulse that produces a selective escape depending 
on the quantum state of the circuit 0, [22[ . The energy 
spectrum versus current bias at $6 = 0.48<I>o and ver- 
sus flux bias at I b = are plotted in FigJS^i and Fig|2]3 
respectively. In the following we will denote v^ m as the 
transition frequency between the states \n a ) and |ra a ), 
with the other mode in the ground state. The first tran- 
sition frequency in Fig|2] is the one of the camelback 
phase qubit, z/qi • With a maximum frequency at zero- 
current bias, the system is at an optimal working point 
with respect to current fluctuations (22[. At higher fre- 
quency, the second transition of the s-mode is observed 
with z/q 2 ~ 2^q X . In the flux biased spectrum the third 
transition Vq 3 is also visible. An additional transition 
is observed at about 14.6 GHz, with a very weak cur- 
rent dependence (Figj2^) but a finite flux dependence 
(FigEb)- It corresponds to the first transition of the a- 
mode, Vq X . The s- mode transition frequencies drop when 
&b/&o approaches 0.7 which is consistent with the crit- 
ical flux $ c /$ = 1/2 + L/(2ttLj) = 0.717 for which 
ujp — »• and x™ in —> tt/2. On the contrary uo® remains 
finite when <fr b ->> <£ c with uo^ ->> y/(4Ej /m)(Lj/L). One 
also observes a large level anti-crossing of about 700 MHz 
between the two transitions v§ 2 an d ^01 • Additionally no 
level anti-crossing is measurable between Vq 3 and z/qi • 

The experimental energy spectrum in Fig J2] can be per- 
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FIG. 2: Energy spectrum. Escape probability P esc versus 
frequency as a function of current bias (a) and flux bias (b) 
measured at <E>b = 0.48<E>o and h = respectively. P esc is 
enhanced when the frequency matches a resonant transition 
of the circuit. The microwave amplitude was tuned to keep 
the resonance peak amplitude at 10%. Dark and bright blue 
scale correspond to high and small P esc - The red dashed lines 
are the transition frequencies deduced from the spectrum of 
the full hamiltonian with C — 510 fF (see text). The green 
diamond is the initial working point for the measurement of 
coherent free oscillations between the two modes, presented 
in FigO and the green dotted square is the area where these 
oscillations take place. 



fectly fitted by deriving the spectrum of the full Hamil- 
tonian described above. As all the other parameters of 
the device can be extracted from switching-current mea- 
surements (22[, the only free parameter is the junction 
capacitance. From the fit we obtain C = 510 fF, which is 
consistent with the 10 fim 2 junction area. Numerical cal- 
culations are used to describe the energy dependence in 
Figl2^i in order to take into account the current bias and 
the inductance asymmetry as well as additional coupling 
terms. The model describes well the level anti-crossing 
between the quantum states |2 s ,0 a ) and |0 s ,l a ), which 
is given by the nonlinear coupling term hg2\X 2 s x a . The 
coupling strength #2i/27r strongly depends on the work- 
ing point. It is predicted to change from zero at $5 = 0, 
to 700 MHz at the anti-crossing, up to about 1200 MHz 
at $6 w 0.65 close to the critical line. The coupling term 
x 2 x a mainly couples the states |0 S , l a ) and |2 S , a ). Start- 
ing from these uncoupled states at <£& = 0, they become 
a maximally entangled state at the resonance condition, 
^01 = ^02- Close to $5 = <1> C the states are still en- 
tangled because #21 diverges. In this device the second 
nonlinear coupling term, x 2 x 2 , is one order of magnitude 
smaller with #22/2^ ~ 50 MHz. Nevertheless the shift 
of transition frequencies that it induces must be taken 
into account to fit the experimental spectrum. Finally 
the absence of x^x a in the coupling term explains why 
no level anti-crossing between ^q 3 and v§ x is observed. 

From spectroscopic measurements, we obtain mini- 
mum resonance widths as narrow as 4 MHz and 3.5 MHz 



respectively. The cor- 
responding dephasing times are estimated to be T| ~ 
160 ns and T% & 180 ns. At the working point (/& = 
0,^6 = Q.37^0) Rabi-like oscillations are performed on 
the s- and a-mode by applying microwave power at the 
resonance frequencies Vq X and z/qi respectively Rabi- 
decay time is measured in the two-level limit with about 
170 ns and 50 ns for s- and a-mode. Relaxation times 
T( = 200 ns and Tf = 74 ns are extracted from the ex- 
ponential population decay of the excited levels |l s ,0 a ) 
and |0 s ,l a ). The measured coherence times of the a- 
mode are much smaller than expected from the minimum 
linewidth. The origin of this additional decoherence will 
be discussed later. Nevertheless, the coherence time is 
sufficiently large for independent coherent control of each 
mode. 

One of the opportunities given by the rich spectrum of 
this two DoF artificial atom is the observation of a coher- 
ent frequency conversion process using the x 2 x a coupling 
of Eqj2j The pulse sequence, similar to other states swap- 
ping experiment 29|, 30[, is presented in (see Figj3^). At 
t = 0, the system is prepared in the state |0 s ,l a ), at 
the initial working point $5 = 0.37<£>o (green diamond 
in FigjJJ}). Immediately after, a non-adiabatic flux pulse 
[32| brings the system to the working point defined by 
3>int, close to the degeneracy point (vq 2 ~ ^01 )• Af- 
ter the free evolution of the quantum state during the 
time Ati n t, we measure the escape probility P eS c Figj3}} 
presents P esc as a function of Ati nt for $ int = 0. 515<£>o. 
The observed oscillations have a 815Mi^z-characteristic 
frequency (inset of Figj3}}) that matches precisely the 
theoritical frequency splitting at this flux bias (red ar- 
row). In Fig|3t, we present these oscillations as func- 
tion of ^int- Their frequency varies with <£>i n £, showing 
a typical "chevron" pattern. In the inset of Figj3h, the 
oscillation frequency versus &i nt is compared to theoret- 
ical predictions. The good agreement between theory 
and experiment is a striking confirmation of the observa- 
tion of swapping between the quantum states |0 S , l a ) and 
|2 s ,0 a ). Instead of the well known linear coupling x s x a 
between two oscillators, which corresponds to a coher- 
ent exchange of single excitations between the two sys- 
tems, here the coupling x 2 x a is non-linear and produces 
a coherent exchange of a single excitation of the a-mode 
with a double excitation of the s- mode, i.e. a coherent 
frequency conversion. Starting from the state |0 s ,l a ), 
an excitation pair |2 s ,0 a ) is then deterministically pro- 
duced in about a single nanosecond at the degeneracy 
point $ int = 0.537^0- 

We now discuss the coherent properties and measure- 
ment contrast in our device. The unexpectedly short 
coherence time of the a-mode can be explained by the 
coupling to spurious two- level systems (TLS) [3l|. With 
a junction area of 10/im 2 our device suffers from a large 
TLS density of about 12 TLS/GHz (barely visible in 
Figl2|). Therefore it is very difficult to operate the a- 
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FIG. 3: Free coherent oscillations between states 
|0 s ,l a ) and |2 s ,0 a ) produced by a nonlinear coupling. 

(a) Schematic pulse sequence. The energy diagram, with- 
out coupling in blue/red and with coupling in black, is rep- 
resented for both the preparation and interaction steps, (b) 
Escape probability P esc versus interaction time At int- The 
inset presents the Fourier transform of these oscillations with 
a clear peak at 815 MHz. The red arrow indicates the theo- 
retically expected frequency (c) P esc versus interaction time 
At int for different interaction flux &i n t close to the resonance 
condition between isq 2 and z/qi • For clarity the data is numer- 
ically processed using 200 MHz high-pass filter (dashed line in 
inset of (b)). Inset : oscillation frequency as function of flux. 
The dashed red line shows the theoretical predictions. 



ing the system at $6 = will lift this limitation since 
it is an optimal point with respect to flux noise. The 
small oscillation amplitude in Fig |3] has two additionnals 
origins. First the duration t n of the 7r-pulse applied for 
preparation of the state |0 s ,l a ) has to fulfill the condi- 
tion t~ l < v^ 2 — z/q-l to avoid multilevel dynamics. Since 
v \2 ~ v oi = 18 MHz in the present experiment, we fixed 
tj[ = 80 ns. However t n is of the order the relaxation 
time of the a-mode which implies a strong reduction of 
the |0 s ,l a ) occupancy after the 7r-pulse, to about 40%. 
The second origin of the small contrast raises the ques- 
tion of the readout contrast between the states |0 s ,l a ) 
and |2 s ,0 a ). In absence of coupling the escape probabil- 
ity of the state |0 S , n a ) should not be sensitive to n a . On 
the contrary, the escape probability of the state |n s ,0 a ) 
is very sensitive to n s . This difference should lead to a 
strong contrast between the states |0 s ,l a ) and |2 s ,0 a ). 
However, since the coupling strength is very large close 
to the critical lines, |0 s ,l a ) and |2 s ,0 a ) are still entan- 
gled when the escape occurs. This leads to an additionnal 
reduction of the contrast between the two states that pre- 
vents the extraction of the theoriticaly high efficiency of 
the frequency conversion [27| . 



In conclusion, we have designed and studied a new type 
of superconducting artificial atom with two internal DoF. 
The spectrum of this device is well described by consid- 
ering two anharmonic oscillators coupled via nonlinear 
coupling terms. Coherent manipulation of the two DoF 
is demonstrated and the strong nonlinear coupling al- 
lows the observation of a coherent frequency conversion 
between the two internal DoF. If inserted in a microwave 
cavity this process should enable parametric amplifica- 
tion or generation of correlated microwave photons. Fi- 
nally by reducing the critical current of the junctions, 
keeping the ratio Lj/L constant, one can increase the 
strength of the coupling term ^22/27r up to 500 MHz 
providing a strong dispersive frequency shift proportional 
to the population of each mode. This should allow the 
realization of C-NOT quantum gate inside the two DoF 
of this artificial atom or enable an efficient readout of 
the camelback phase qubit state by probing the transi- 
tion frequency of the a-mode. Considering the long list of 
unique and interesting features, this new artificial atom 
with two DoF and a V-shaped energy level structure is 
a valuable addition to the growing set of building blocks 
for superconducting quantum electronics. 



mode in a frequency window free of TLS since is only 
slightly flux dependent. However this is not a real is- 
sue as it can be solved easily by reducing the junction 
area [33]. The minimum linewidth of both a-mode and 
5-mode, and therefore their coherence times, are limited 
in our experiment by low frequency flux noise. Operat- 
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